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Generalized parton distributions for qq pions 

B. C. Tiburzi and G. A. Miller 
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(Dated: February 1, 2008) 

We use a two-body, light-cone pion wave lunction (which vanishes when one of the constituents 
carries all of the plus-momentum) to construct a double distribution. We show that the resulting 
generalized parton distribution has incorrect behavior at the crossover between kinematic regions. 
Furthermore the positivity constraints are not respected. On the other hand, the field-theoretic, 
two-body double distribution gives generalized parton distributions that have physically correct 
crossover behavior and respect the positivity constraints independent of the model wave function 
used. However, the sum rule and polynomiality conditions cannot be satisfied with a two-body wave 
function alone since Lorentz symmetry is not maintained. 
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I. INTRODUCTION 

In recent years generalized parton distributions 
(GPDs) |l], |2|, H) have generated a considerable amount of 
attention. These distributions stem from hadronic ma- 
trix elements that are both non-diagonal with respect 
to hadron states and involve quark and gluon operators 
separated by a light-like distance. Thus physics of both 
inclusive (parton distributions, e.g.) and exclusive (form 
factors, e.g.) reactions is contained in the GPDs. At the 
leading-twist level, these new structure functions describe 
the soft physics of a variety of hard exclusive processes 
(see the reviews H). 

Recently two-body light-front wave function models of 
the pion have been used to obtain GPDs Q based on dou- 
ble distributions (DDs) @, 0. Such wave functions, based 
on solving the Weinberg equation |§] and carrying the 
total charge of the pion vanish when either the quark or 
anti-quark carries all of the plus-momentum of the pion. 
The purpose of this note is to show that the use of such a 
model leads to GPDs with undesirable features. There- 
fore we believe that constructing simple models with the 
relevant physics that satisfies the reduction relations for 
GPDs remains a considerable challenge. 

We proceed using the following logic. The beginning, 
in section O, discusses the pion model GPD introduced 
in H by reviewing the power-law wave function used and 
showing how the DD was obtained from the expression 
for the form factor. We also indicate the reduction rela- 



tions required of the GPDs. In the following section (III), 
the physics at the crossover (X — () between kinematic 
regions is reviewed in terms of the Fock space representa- 
tion. Internal consistency demands a model with a van- 
ishing distribution function (at zero plus-momentum) to 
have vanishing GPDs at the crossover. The proposed 
model is shown not to satisfy this restriction. Further- 
more, the realistically modified model of Ref. (5) docs 
not match up to the physical intuition afforded by the 
Fock space representation which relates small x behav- 
ior of q(x) to 
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GPDs at the crossover. Section 
vestigates whether the models respect the positivity con- 
straints. The realistically augmented model satisfies pos- 



itivity, while its underlying wave-function-based model 
violates the constraint for all X > £. The Fock-space 
representation is then used to obtain DDs by using the 
valence wave function consistently in the relevant ma- 
trix element (section |v|). These DDs are shown to violate 
Lorentz symmetry by retaining ^-dependence. The two- 
body model GPD obtained from these DDs vanishes at 
the crossover and satisfies the positivity constraints. As 
a result of the ^-dependence, however, moments of these 
GPDs do not satisfy the polynomiality conditions. In 
particular the sum rule cannot be satisfied in the valence 
sector alone. Additionally in the Appendix, we scrutinize 
the covariant scalar model used in the Appendix of to 
justify the construction of the DD. Parametric differen- 
tiation of a covariant field-theoretic amplitude is shown 
not to preserve positivity. A brief summary is presented 
in section 



VI 



II. PION MODEL WAVE FUNCTION 

Here we review the basics of the model proposed by 
H starting from the two-body wave function and ending 
at the double distribution. Motivated by a numerical 
solution to an effective two-body equation for the pion 
M, the following wave function is adopted 



N 



^{l^x)[a{x)+b{ X )^Y 
where a(x) and b{x) are the dimensionless functions 



a(x) 



1 



'■(x- 



1)2 



x(l — x) 



b(x) 



4x(l - x) 1 



(1) 



(2) 



with s a dimensionless parameter of the model and N 
the normalization constant. One then determines the 
normalization N from the condition 



J dxdk 1 - ^(x,^)] 2 = 1, 



(3) 



which makes us think of tp as a wave function that in some 
way effectively represents contributions from higher Fock 
components. 
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Having modeled the physics in the lowest (gosymmet- 
ric) Fock component, the Drell-Yan formula || is used 
for the form factor 

F ( t ) = /^.(^ +(1 -*)A^(^), (4, 

where t = — A ±2 . Using the explicit structure of the 
wave function Eq. (|l]), the form factor can be manipu- 
lated into the form 



In this form, the sum rule for the form factor 



F(t) = [ dx [ 
Jo Jo 



dy F(x,y;t), 



(5) 



with the DD 



F(x,y;t) = 0(l-x-y) 



6y{l-x-y)/(l-x) 3 



(6) 



Here A is just a constant related to the normalization by 

_3_ s 2 1 

47r m' 2 TT 2 



A = 1:^4. This DD has the proper support and is 
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Miinchen symmetric |1CJ]: y; t) = 1 — x — y;t). 
In the forward limit, the y-dependence factorizes 



F(x, y; 0) = 0(1 - x - y)h(x, y)q(x) 
into the normalized profile function h(x, y) 

6y(l -x-y) 



h(x,y) 



(l-x)2 



satisfying 



l-x 



dy h(x,y) = 1. 



(7) 



(8) 



(9) 



The profile function h(x, y) is identical to the asymptotic 
quark distribution amplitude model suggested in . The 
properties Eqs. (Q) and (^|) are essential so that the total 
charge F(Q) looks like an integral of the quark distribu- 
tion function q(x), namely 



q(x) = J dk 1 - \ip(x,k^)\ 



Aa(x) 



(10) 



Notice that q(Q) = 0. 

The ingenuity of DDs comes about when we construct 
the GPD H via 



Tc_{X,t) = ( dx I dy F( 
Jo Jo 



x,y;t)S(X -x-Cy). (11) 



f dX F c (X,t) = F(t) (12) 
Jo 



and the polynomiality constraints on moments of T are 
immediately satisfied. 

In summary, the above model satisfies the spectral 
properties and reduction relations required of GPDs. 
Moreover, unlike previously used Gaussian Ansdtze for 
C = GPDs, there is non-trivial interplay between X, £ 
and t dependence. 

III. AT THE CROSSOVER 



The first inconsistency of the qq pion model comes 
from behavior at the crossover, i.e. when X = £. Using 
the Fock space representation of deeply virtual Compton 
scattering |ll], [TJ| , we first derive a result that focuses on 
the small-x behavior of the model wave function. 

Since q(0) — above, let us consider models such that 
the quark distribution function q(x) vanishes at x = 0. 
This quark distribution function, however, is a model 
of all the higher Fock components. Although we can- 
not generally disentangle these components from a given 
model for ip, we implicitly have the restriction 



n j "'{"I \i 



(13) 



where the sum runs over all j-quarks in each n-body Fock 
state and the integration measure is given by 



_ (14) 

But since each term's integrand in Eq. dTJf) is positive 
definite, we must have 



ip n (. ..,Xj=0,...)=0 



(15) 



whenever j corresponds to a quark. 

This strong result (which applies only to models with 
q(Q) — 0) has a direct implication for GPDs at the 
crossover. Approaching the crossover from above we have 



^(Ct)cx^(i-c) 1 -^^ f ^sic-xM'ix'iX^WMzirt,*), 



(16) 



with the primed variables given by 



Thus since the struck quark rebounds with zero plus- 



4 = 



i-C 



= kf - A 1 - (17) 
k^-kl+^A 1 . (18) 
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FIG. 1: Value of the generalized parton distribution at the 
crossover: t). Here we choose the various values for £ 

and plot the crossover value as a function of ■=£ for the model 
parameter s = 0.95. 



momentum, Eq. ( |l5| ) forces T^{C,,t) to vanish. 1 Con- 
sistency requires any effective modeling of the higher 
Fock components with a vanishing quark distribution at 
zero plus-momentum to mandate GPDs to vanish at the 
crossover. Here we have considered the GPD for a scalar 
particle. For a spin one-half bound state, the correct 
generalization of the above result is that the imaginary 
part of the deeply virtual Compton amplitude M van- 
ishes. This is because the generalized parton distribu- 
tions H(X, C, t) and E(X, £, t) enter as the linear combi- 
nation 



1 2 



( 2 



4(1 -§yw 



(19) 

which has a light-cone Fock space decomposition similar 
to Eq. (fl6|); the only differences are labels for identical 
initial and final bound-state spin. 

The pion model of section || has a quark distribution 
function ( [l0| ) which vanishes at zero plus-momentum. 
Thus the model effectively takes each Fock component to 
vanish at x — 0, and hence we would expect this model's 
GPD (X, t) to vanish at the crossover. Using Eq. jll| ) 
for the model DD (^|), we find it is not the case (as one 
could already see from Figure 4 in [El) 



*c(C,f) 



dyF(((l-y),y;t) 



(20) 



In Figure we show the value at the crossover as a 
function of for a few values of (. Despite the vanishing 
of the quark distribution function at zero x, the model 
GPD of [pj does not vanish at the crossover. Said an- 
other way, the non- vanishing behavior at the crossover is 



1 There is the possibility of the Gibb 's phenomenon, for which the 
sum of infinitely many zeros in Eq. ( |16[ ) need not vanish. This 
is not expected; Fock components are not independent. 



el parameter is again 
we expect the original 



FIG. 2: Comparison of the model GPD Eq. QllJ) (MODEL) 
with the realistic, factorized form Eq. (p|) (FACTORIZED) 
for fixed £ = 0.4 and t = -m 2 . The mo. 
chosen to be s = 0.95. Physically 
model to vanish at X — ( while the realistic version should 
be greatly enhanced there. 



indicative of higher Fock space contributions which are 
not present in the quark distribution function, i.e. the 
constituent quarks in Eq. (Q) have no sub-structure. The 
model GPDs appear to have more physics than the un- 
derlying wave function. 

Acknowledging the limits of the model quark distribu- 
tion function (10), the authors || suggest implanting a 
more realistic distribution 



q R {x) 



3 1 

4 "T 



(21) 



and modify the DD accordingly by using a factorization 
assumption 



F R {x,y-t) = F(x,y;t) 



q R {x) 
q(x) ' 



(22) 



Modifying the quark distribution function according to 
Eq. ( pl| ) must change the wave function. Consider- 
ing the singular behavior, the effective wave function 

ipeff — must be dramatically enhanced at small-x. 

Thus one consequence we expect is a sizable enhance- 
ment to the GPD at the crossover due to the form of Eq. 
(|l6|). Comparing the original GPD (|IT| ) to the realistic 
one calculated from the DD in Eq. ( ]22[ ) in Figure ||, we 
notice that the integral of Eq. ( P0|) leads to only a small 
enhancement. Thus the factorized assumption does not 
support an overlap interpretation in terms of the effective 
wave function. 

One could ask why bother looking at the crossover in 
the first place. Experimentally [jl3| the DVCS amplitude 
is accessed from the dominant Bethe-Heitler contamina- 
tion through beam asymmetries. At first the interfer- 
ing QED radiative processes were seen as a hindrance 
to unmasking the DVCS contribution to the total cross- 
section. Now however, they seem more of a virtue, al- 
lowing access to both the real and imaginary parts of the 
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virtual Compton amplitude. The beam-spin asymmetry 
for example, allows the imaginary part of the am- 
plitude to be measured. At leading twist, the imaginary 
part involves the GPDs evaluated at the crossover X = (. 
From the point of view of modeling, it would be nice to 
understand measurements of the beam-spin asymmetry. 
The above model, however has been constructed || so 
that it should yield a vanishing value at the crossover. 
Any non-zero result of a qq model, is simply an artifact 
of Eq. (||). As we shall show below in Sect. [v| the inte- 
grand of that expression is not the DD because it is not 
obtained from is definition as a specific matrix clement 
of field operators. A better determined qq wave function 
will not lead to a better understanding of GPDs at the 
crossover. 

On the other hand, if one is concerned with the real 
part of the virtual Compton amplitude (which enters in 
the charge asymmetry |H| ) , the GPDs appear in a (prin- 
ciple value) integral weighted by the hard scattering am- 
plitude. The crossover effects are still important because 
a factor of is present and the GPDs are likely to 
have discontinuous derivatives. The only sophisticated 
model study on the market |]l6f (which arguably does 
QCD in the large N c limit) shows considerable contribu- 
tion (~ 60%) to the cross-section from near the crossover. 
Lastly one could hope to measure the GPDs directly via 
double DVCS fil. 



IV. POSITIVITY CONSTRAINTS 

The result proved in the light-cone Fock representa- 
tion in section 



[II concerning behavior at the crossover 



stemming from the small- a; limit of the quark distri- 
bution also follows from the positivity constraints for 
GPDs. Originally these constraints appeared in |?], [l^] 
and were derived from the positivity of the density ma- 
trix by restricting the final-state parton to have posi- 
tive plus-momentum (and ignoring the contribution from 
E(X,(,t) for the spin-i case). Although the matrix 
elements involved for GPDs are off diagonal, they are 
still restricted by positivity and their diagonal elements. 
Correcting the constraints for the presence of the E- 
distribution was first done in |fj"2f . By considering the 
positivity of the norm on Hilbcrt space, stricter con- 
straints for the spin-i distributions H and E have re- 
cently appeared as well as constraints for the full set of 
twist-2 GPDs @. 

For the scalar pion case there is of course no contribu- 
tion from the non-existent E'-distribution and hence the 
original bounds are actually correct. When X > ( the 
positivity constraint for the spin-0 pion reads 



R(X,Q = (l-C/2) 



\T c (X,0)\ 



(23) 



Of course the result holds for finite —t, however since the 
function T decreases with —t, so that (|23|) is the tightest 
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FIG. 3: Comparison of the model GPD Eq. (|llj) (MODEL) 
with the realistic, factorized form Eq. @ (FACLORIZED) 
for fixed ( — OA a,t t — 0. The model parameter is again cho- 
sen to be s — 0.95. Here we plot the ratio R(X, Q appearing 
in equation ( |23| ) as a function of X > £. Positivity constrains 
this ratio to be less than one. 



constraint. Notice the bound at X = ( corresponds to 



the result derived in section [II provided of course the 
hypothesis is met: q(0) = 0. The GPD derived from Eq. 
(p) violates the constraint for all X > £ as depicted in 
Figure || for ( — 0.4, while the augmented model fl2^ ) 
satisfies the constraint. 

In the realistic quark distribution model ( pl| ) the pos- 
itivity constraint is quite lenient as q R (0) — oo allows 
for no bound on the GPD at X — £ in equation (|2^). 
Consequently near X = ( the GPD is only very loosely 
constrained. In general, a small- a; factorized fix-up of 
q(x) will not guarantee that the positivity constraint will 
be met. 



V. DOUBLE DISTRIBUTIONS IN FOCK SPACE 

As we have seen the pion model proposed in || violates 
the positivity constraint near the crossover. The realistic 
fix-up, while satisfying positivity, does not agree with the 
intuition provided by the light-cone Fock space represen- 
tation which, e.g., connects small- a; behavior to GPDs at 
the crossover. Our solution to this dilemma is to define 
the double distribution consistently using the same Fock 
space representation implicitly used in Eqs. (|l(]) and (Q). 
To achieve this, we truncate the Fock space at the valence 
sector as we did previously |2(J . 

In the two-body truncated Fock space the pion bound 
state of momentum P appears as 



\ir(P) 



dxdk 



^J1QtT 3 x(1 — x) 



V(l,2) &t(i )d t( 2 )|o>, 



(24) 

where 1 corresponds to the momentum xP + , k- 1 , while 2 
represents (1— x)P + , P -1 — k ± . Without loss of generality, 
we have chosen particle 1 to be the quark and 2 the anti- 
quark. One can use the theory's Galilean invariance to 
write ip (1,2) = ip(x, k- 1 — xP ), which is the same as 
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— XjxP 1 - — k -1 ) by qq symmetry. We omit helicity 
labels to make our analysis parallel to that of ||. The 
bound state is normalized according to 

( tt(P') I tt(P) ) = 16ir 3 P + S(P ,+ - P+)(5(P' ± - P x ), 

which in turn implies the following normalization for the 
valence wave function 



J dxdk 1 - |^(cc,k J 



1. 



As above, we can think of Eq. (g) as the consequence 
of truncating the Fock space to only the valence sector, 
i.e. ip must be an effective two-body wave function. The 
constituent quarks in our bound state, however, have no 
sub-structure nor should we gain Lorentz invariance from 
Eq. (|). 



The good component of the quark field tp + = V+ip 
(with P + = 7~7 + /2) has a mode expansion 



X ) 



dk+dk J -9(k+) 

167T 3 fc+ 



M)V+u x (k+M)e- %k+z ~ 



d t (fc+,k x )P +UA (fc+,k ± )e lfc+J 



(26) 



where the mode operators satisfy the anti-commutation Additionally the matrix element definition of the form 



relation 

\b{k+, k ± ), &t(p+, p-L) J = 16^ 3 fc+<5(fc+ -p+)<5(k ± -p x ) 

and similarly for d and (V . 

These equations are all consistent with the starting 
point of the pion model in section ||. For instance: the 
quark distribution function is identical to Eq. ( |To| ) 

q{x) = lj ^- e ™ p+ *-{n(P) |tf(0) 7 + V(OI tt(P)>. 

(28) 



factor lines up with the Drell-Yan formula 



F(t) 



2P+ 



(tt(P') |^(0)7+^(0)| tt(P)), (29) 



which is given by Eq. (0) with P' = P + A and t = A 2 = 
-A ±2 , i.e. A+ = 0. 



J 



A. Derivation of the DD 



The DD is defined from the non-diagonal matrix element of bilocal field operators 01- To derive the DD 
consistently, we should respect this underlying field theoretic construction: 



2P+7I - C/2) MP/) IV>(°h + ^(OU(^)) = y o J q 9(l-x-y)F(x,y;t) 



e -i{xP + -y& + )z- _ e i{xP + + (l-y)A^ 



Using the truncated Fock space expansion of the bound state (E4h, we can express the matrix element as 



dxdy. 
(30) 



(tt(P + A) | ^(0)7+^(0 I k(P)) = 2 / dk+dk 1 - 9(k + + A+)6»(fc+) 

x ip*(Y,\s x +(l-Y)A ± )ip(X,lc 1 -) 



-ik + z 



, (31) 



where we have used the temporary replacements: X = k + /P + and Y = -pr+^r- To find the DD, we merely 
need to cast the above into the form dictated by Eq. (|30|). Thus it remains to introduce x,y. Before doing so 
explicitly, it is wise to make a few comments. Given the z~ dependence of Eq. (|3l|), we must introduce x,y only via 
k + = xP + — yA + . In turn, however, the resulting F(x, y; t) must depend on A+ = — (P + . This is obvious looking at 
the two-body partonic cartoon shown in Figure [I] for the DD. In fact, one can see that any iV-body Fock component 
contribution to the DD will retain £ dependence. This should not be surprising: the light-cone wave functions are 
not covariant objects, ft is only when summing over all Fock component overlaps that the £ dependence in the DD 
cancels. 
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FIG. 4: Partonic diagram for valence two-body double distribution. 



To cast Eq. ( pl| ) into a form from which we can read off the DD, we must use an explicit form for the wave function. 
Using the model wave function (|l|) above, we can write 

m z N £ J 



dk A 



[a(X) + b(X) k^ 2 ] 2 [a(Y) + b(Y) ( +(1 - Y) A^ 2 )] 2 ' 



(32) 



Notice that now k^ and A have been scaled by m to be dimensionless. Exponentiation of the denominators via 
A~ 2 = ae~ aA da allows for a Gaussian integral over k 1 - to be performed. The result is 

°° r °° «i«2 J /v , , v s a ia2 b(X)b(Y)(l-Y) 2 A ±2 ^ 



I = n / da 1 da 2 — u —— exp < - a>ia(X) - a 2 a(Y) — — — — \. (33) 

Jo Jo aib(X) + a 2 b(Y) aib(X) + a 2 b(Y) j 

Since b > 0, we define the new variables X,(3: A = a\b(X) + a 2 b(Y) and a.\b(X) = A/3. Performing the resultant 
integral over A yields 



[b(X)b(Y)f J 



/3(l - p)d(3 



'b(X) 



(34) 



(y. Under this change of variables 9(f3)6(l — (3) 



Next we introduce y and a by y — (1 ~ a 2 )P and X = a 2 

9(y)9(l — y)0(l — a 2 )9(l — a 2 — y) and further d(idk + = 2P + |a| ^" a 7^ dyda. The resulting integral is even in a 

and hence we introduce x — a 2 so that J , da\a\ . . . — J Q dx Noticing that 9(1 ~ a 2 ) = 9(x)9(\ — x), we can write 

out the DD from Eq. M 



F(x,y;t) 



6s 2 (l - - X)9(X - C) 



y(l-x-y)/(l-x) 4 



4A(1 - C,/2)y/X(\ - X)Y(l - Y) b 2 (X)b 2 (Y) 

le labels are X = 
clearly different from Eq. (^J). 



y a W i l-x-y a ( Y ) I y(i-^-y) 



(i -r) 2 a j 



(35) 



where now the labels are X = x + £y and Y" = x ■ This DD based on the field-operator matrix element (|3(]) is 



B. Properties 

The use of a qq wave function is a violation of Lorentz 
symmetry which has drastic consequences for the com- 
puted DD of Eq. d35|). The violation of Lorentz sym- 
metry and the Miinchen symmetry is manifest in the in- 
dependence of F in Eq. (|3^). This C-dependence spoils 
the sum rule for the electromagnetic form factor. Addi- 
tionally the DD is restricted by C in x— y space as a result 
of the ^-functions that maintain the light-cone spectrum 
condition. One way to obtain the full support proper- 



ties of the DD is to use higher Fock components. Alter- 
natively, one could imagine constructing a C-dependent 
function for the excluded region of phase space which has 
the properties necessary to cancel the ^-dependence of F 
above. In essence, this would be an attempt at solving for 
the (light-cone) instantaneous kernel for -0 and utilizing 
crossing to access all kinematical regions One can 

do so in the case of an instantaneous interaction because 
there are no true higher Fock states. 

Taking the ( — > limit of Eq. ( |35| ) results in the 
pion model DD above (^). Thus we additionally have 
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the reduction relation at t — (Q). It is not surprising 
that the A + = limit corresponds to the model of ||. 
They extract the DD from the Drell-Yan formula which 
translates to calculating the matrix element in ([50|) using 
the A + = frame. It seems reasonable that calculation 
of the Lorentz invariant DD can be done in the A + = 



J 



frame; however, the result does not satisfy positivity. 

Using Eq. (p0|), we can show that the GPD resulting 
from Eq. (|3^) vanishes at the crossover. This follows 
simply from the Y — > behavior of F. Working with 
wave functions directly illuminates this. From Eq. d3l|), 
we can write down the GPD 



6{X)6{1 - X)0{X - CW(Y, k x +(1 - Y) A ± )ij}(X, k x ) 
X)6(X - c + 1)6(( ~ x)r ( - X/ {I -C)M +(1 + x/(i - c) a 1 - ) v(C - x, k x ) 



(36) 



X 



0.2 




and the particular form of Eq. (gTJ) leads to a singular 
GPD at the crossover. The effect of qn(x) is to mimic 
contributions from higher Fock states. Continuity of the 
GPD arises from relations between Fock components at 
vanishing plus-momentum. A pcrturbative model con- 
firms this [ p2[ . Constructing simple light-cone wave func- 
tion models for GPDs consistent with the reduction re- 
lations and positivity remains a challenge. 



VI. SUMMARY 



FIG. 5: Plot of the positivity constrained GPD Eq. (|36j) for 
£ = 0.4 and with the model parameter s — 0.95. We choose 
A x = 0.6m 2 so that t will be essentially the same as in 
Figure 



From the above expression, we directly see the GPD's 
vanishing at the crossover results from the small- a; be- 
havior of the model wave function. Furthermore, as is 
generally true |L2| , direct application of the Schwartz in- 
equality yields the positivity constraint (23]). In Figure 
|, we plot the GPD Eq. ( pf) w hich can also be obtained 
from (fv|) by means of E_q. (|11|) . There is a big difference 
between Figures [| and Q 

Lastly one could try to incorporate a realistic quark 
distribution by using a factorized Ansatz for the double 
distribution. Implementing Eq. (|2^) for the DD above 
(|35|), however, does not affect the vanishing of the GPD 
at the crossover. This is again contrary to our observa- 
tion from Eq. ( |l6| ) and hence the factorized assumption 
must not change the small- a; behavior of the wave func- 
tion. A more physical Ansatz is effected by directly alter- 
ing the behavior at the end-point: iji — » ip(x, k x )y^^. 
This alteration maintains the positivity constraint re- 
gardless of the functional form of q^(x). The drawbacks 
of this form are clear: resulting GPDs are discontinuous 



The DD based model of GPDs proposed in [|| is ap- 
pealingly ingenious, but has been shown to have undesir- 
able properties. Behavior at the crossover is inconsistent 
with intuition. Positivity is violated because additional 
time-orderings involving non-two-body contributions are 
omitted. Alternately we can view violation of the positiv- 
ity constraint as resulting from parametric differentiation 
of a covariant field-theoretic amplitude. As shown in the 
Appendix, such differentiation does not preserve positiv- 
ity. On the other hand, the light-front Fock representa- 
tion of GPDs is intuitively physical, satisfies positivity, 
but manifestly breaks Lorentz symmetry. 

We present these problems for two reasons. Firstly, 
to treat other theoretical attempts in this direction with 
caution. For example, one could imagine applying the 
procedure in Q to a better determined (or better tuned) 
valence wave function. We have shown that this will 
not result in a more accurate determination of the GPD. 
The procedure is not internally consistent (i.e. the pos- 
itivity constraints will be violated) and one must ulti- 
mately augment the wave function with a realistic quark 
distribution using an untested factorization Ansatz. This 
augmentation is an attempt to alter the small- a; behav- 
ior of the wave function. However, the resulting GPDs 
do not match up with the physical intuition provided by 
the light-cone Fock representation at the crossover which 
closely relates the low-cc quark distribution functions to 
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GPDs at the crossover. Secondly we believe this cau- 
tionary note could be important to the experimentalists 
hoping to interpret future data. The reverse of this pro- 
cedure seems to suggest one can learn about the effective 
valence wave function from experimental access to GPDs. 
There is no such correspondence: the model GPD from 
Eq. (^|) appears to contain more physics than the valence 
wave function from which it was derived. 

Lastly we remark that the proposed model |j| does 
provide an alternative to Gaussian based forms — allowing 
for the correct ^-dependence. While the model is useful 
as an analytical result, one cannot take the X and £ 
dependence seriously. 



Acknowledgments 

We are indebted to A. V. Radyushkin for discus- 
sion and correspondence. This work was funded by 
the U. S. Department of Energy, grant: DE-FG03 — 
97ER41014. 



APPENDIX A: DOUBLE DISTRIBUTIONS IN A 
SCALAR MODEL 

Above we have seen that the wave function based DD 
(^) does not satisfy the positivity constraint which cre- 
ates doubt about the uniqueness of its construction. In 
their Appendix, the authors Jpj, however, do not treat the 
wave function as in section |V|. Instead, they derive the 
DD in a covariant fashion from the handbag diagram with 
point-like couplings. By taking parametric derivatives 
with respect to the initial- and final-state mass-squared, 
Eq. (||) is derived for this model. This is the justification 
presented that the procedure outlined in section || leads 
to the correct DD. Below we investigate positivity in this 
model and show that violation of the constraint comes 
about because of the parametric differentiation. 

Using the one-loop handbag diagram with scalar cur- 
rents in the Bjorken limit, the vector-current DD is de- 
rived n 

F v (x,y;t,pl,pl) = xi-y{l-x- y)t 

1 -l 

-x[(l-x-y)pl + ypl] + m 2 j , (Al) 

where p\ and p 2 , are the initial- and final-state pion four- 
momentum squared, respectively. Here they are treated 
as free parameters. Evaluating at the pion mass M 2 , we 
arrive at the DD 



m 2 — x(l — x)M 2 — y(l — x — y)t 
with a corresponding quark distribution function 



x(l — x) 



m? — x(l — x)M 2 



(A2) 



(A3) 



To derive the DD encountered in section ||, one applies 
{p\d / dp\){pld / dpi) to F v in Eq. (|Al|). Evaluating at the 
pion mass, they find |5| 



F ( v\x,y;t) 



2y(l-x~y)x 3 M 4 



[m 2 — x(l — x)M 2 — y(l — x — y)t] 3 ' 

(A4) 

with a corresponding quark distribution 



(2)/ x M 



x 3 (l — x) 3 



3 [m 2 - x(l - x)M 2 ] 



213 ' 



(A5) 



Eq. (A4) can be manipulated into the form of Eq. 



(up to the overall normalization) using the replacement 



Y = -A4 2 /4: 



Thus one has a covariant way to 



derive power-law DDs. 

To investigate positivity, we must first derive the 
proper constraint. Neither models satisfy Eq. (p3|), which 
is surprising for Fy^ since with point-like quark-pion cou- 
pling, there are no excluded non-two-body diagrams in 
the calculation. Eq. ( p3| ) holds for a composite scalar 
of spin-i constituents. The constraint for scalar con- 
stituents is different as we shall now see. 

For a composite scalar of scalars, a matrix element 
definition 2 of (X, t) is 



f c {X,t) 



(2X - ()P+ f dz 



2-C 



2vr 



e iXP+z '(P'\q(0)q(z-)\P) 

because it properly reduces to the quark distribution 
function 

T (X,0) = XP+J ^-e* XP+z -(P\q(0)q(z~)\P) = q(X) 
and integrates to the form factor 



J dX T c {X,t) = 



(P'|g(O)i0+g(O)|P) 
^+(2-0 



F(t). (A7) 



Similar to [|18| , we can derive the relevant constraint by 



using translation q(z 



q(0)e 



and writ- 



ing the Cauchy-Schwartz inequality in the following form 

2 



X (P\q(0)\Z 



'p'\m\z) 



V i-c 

x S(XP+ + p+ - P+) > (A8) 



for the complete set of states {!-£)}• Minimization with 
respect to the parameter a yields the constraint 



x(x-o 



2-C 



*c(*,o) I < i- 

(A9) 



We thank M. Diehl for helping to straighten out this definition. 
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FIG. 6: Comparison of the scalar model, DD-based GPDs: 
Eq. (|A§) (F v ) with Eq. (Fy) for fixed < = 0.4 at t = 0. 

The model parameters are chosen to be m = 0.33 GeV and 
M = 0.14 GeV. Here we plot the ratio R'(X, £) appearing in 
equation ( |A9| ) as a function of X > £ . Positivity constrains 
this ratio to be less than one. 



We can use Eq . (jA9|) to test whether GPDs constructed 
from DDs Eqs. (|A2|) and ( A4) satisfy positivity. To con- 
struct their respective GPDs, we use equation IO. In 
Figure [| we compare the two GPDs. Not surprisingly 

Fy satisfies the positivity constraint. As we commented 
above, there are no missing time-orderings in the calcu- 
lation of this GPD since the point-like coupling between 
the two quarks and pion is covariant. 



As we see from the Figure, the parametrically differen- 

(2) 

tiated DD F v ' violates positivity for all AT > £. When 
one treats the quarks as spin-^ particles, it is thus nat- 
urallyexpected that the GPD calculated from the DD 
Eq. @ violates Eq. (|3|). Viewed in this way, however, 
it is more of an accident that the GPD calculated from 
Eq. (E3) satisfies the constraint. 
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